Quantum mechanical picture of the coupling between interlayer electronic excitations 
and infrared active phonons in bilayer cuprate superconductors 
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The formula frequently used to describe the c-axis infrared response of the coupled electron-phonon 
system of bilayer cuprate superconductors and providing inrportant insights into the physics of these 
materials has been originally obtained at the level of the phenomenological multilayer model. Here 
we derive it using diagrammatic perturbation theory. 
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The phase diagram of underdoped high-T c cuprate su- 
perconductors is divided into three regions: the super- 
conducting one below T c , the pseudogap phase in the 
temperature range T c < T < T* , where T* is the on- 
set temperature of the (high energy) pseudogap, and the 
anomalous normal state above T* X& There is an inter- 
esting and important class of physical properties, that 
begin to develop/ whose temperature dependence changes 
its character at an onset temperature much lower than T* 
but at the same time considerably higher than T c . 3 These 
include, e.g., the NMR relaxation rate (1/TiT)^ ampli- 
tude of the transverse plasma mode (also referred to as 
the 400 cm -1 mode) and some phonon structures in the c- 
axis infrared conductivity of YBa2Cu307_,5 (Y-123)— ~—, 
in-plane scattering rate and in-plane conductivity below 
ca 700 cm -1 ^ii£ ; amplitude of the neutron resonance 11 , 
Nernst signal^, diamagnetic responsei 3 -, specific heat^, 
some features of the excitation gap 15 ' 16 , spectral weight 
at the chemical potential^, magnetic field dependence 
of the in-plane resistivity^, THz conductivit y 19 ' 20 , and 
the quasiparticle interference 2 ^. The values of the onset 
temperature range from ca 15 K above T c (THz conduc- 
tivity) to more than 100 K above T c (transverse plasmon 
and the phonon structures, Nernst). Some of the phe- 
nomena (increase of the THz conductivity and quasipar- 
ticle interference) are almost certainly related to super- 
conducting fluctuations. Some other with a much higher 
onset temperature (e.g., the transverse plasmon and the 
phonon structures, Nernst), have also been attributed to 
a precursor superconducting state (PSC) lacking the long 
range phase coherence. This exciting issue, however, is 
fairly controversial and no consensus has been reached 
thus far. There are two prerequisites for an ultimate un- 
derstanding: (a) A detailed knowledge of both doping (p) 
and temperature dependence of the relevant quantities, 
(b) An understanding of the physical contents of these 
quantities, in particular, the relation to superconductiv- 
ity. 

Dubroka et alr^ have studied in detail the tempera- 
ture dependence of the phonon structures in the infrared 
c-axis conductivity of Y-123 for a broad range of p of 
0.03 < p < 0.17 and determined the p-dependence of the 
corresponding onset temperature T° ns . The relation be- 



tween the phonons and the electronic degrees of freedom 
has been understood in terms of the phenomenological, 
essentially classical, multilayer mode l 22 ' 23 , involving the 
local conductivity of the bilayer units (denoted here as 
o"bi) and that of the spacing layers (<7i n t), extende d 24 ' 25 by 
including formulas for the local fields acting on the rele- 
vant ions and a feedback of the phonons. For a schematic 
representation of the electronic part of the model, see 
Fig. 1 (b) of Ref. [26|. For a careful discussion of the role 
of the phonons, see Ref. [25|. The model has allowed the 
authors of Ref. 10 to conclude that at T ons the spectral 
weight of the low-energy component of <7bi (the Drudc 
peak in the analysis of Fig. 1 of Ref. llOl ) begins to in- 
crease, signalling an increase of coherence of the response. 
Several arguments have been advanced suggesting that 
this is due to the formation of a PSC. The importance of 
the topic of PSC and of other interesting fields of appli- 
cation of the multilayer modeler— calls for developing a 
fully quantum mechanical description of the response of 
the coupled electron-phonon system, that would support 
the phenomenological approach. Chaloupka, Bernhard, 
and Munzar have already constructed (see Ref. [26j, CBM 
in the following) a microscopic gauge-invariant theory of 
the c-axis infrared response of bilayer cuprate supercon- 
ductors, such as Y-123, providing both a justification and 
limits of the electronic part of the multilayer model and 
of a related earlier approach 31 . The phonons, however, 
have not yet been incorporated. Here we provide and ex- 
tension of the theory of CBM, where they are included at 
the Green's function level. A similar approach has been 
recently used to explain the electric field dependence of 
the infrared spectra of bilayer graphene 32 . The formula 
for the c-axis conductivity will be derived and it will be 
shown to coincide, for a natural choice of the electron- 
phonon coupling term, with that of the phenomenological 
model. The essential aspects of the theory are formulated 
in the points (i)-(viii). 

For the sake of simplicity we neglect the hopping 
through the spacing layer (t' ± in the notation of CBM, 
Cint = 0), focus, as usual, on the long- wavelength limit, 
q = 0,— and limit ourselves to the case of one infrared 
active phonon. 

(i) Relevant equations of CBM. The (electronic) con- 



ductivity a b i of the bilayer unit is given by 

o-bi(w) = -zweoXbi(w) = — — j-^ , (1) 

l[U) + 10) 

where %bi is the bilayer susceptibility, Kdia-bi the term 
describing the diamagnetic contribution, dbi the dis- 
tance between the closely spaced planes, and Hjj(u>) the 
current-current correlation function. The term Kdia-bi is 
given by 



3-^r Vi±(k)(4 ks c Bks -c^ ks CAks), (2) 



K " 1_ h 2 Na 2 ^' 

K.S 

where TV is the number of surface unit cells per copper- 
oxygen plane, a is the in-plane lattice parameter, the sum 
runs over all values of the in-plane quasimomentum k and 
spin s, t±(k) is the intra-bilayer hopping parameter, and 

c sks' c Bks ( c Aks> c Aks) are the quasiparticle operators 
corresponding to the bonding (antibonding) band. The 
correlator TLjj is given by 



%(<*>) = 



iNa 2 



dt - lut 



<[&(*),&(0)]Wt), (3) 



where j^ is the operator of the intrabilayer paramagnetic 
current density, 



Jbl 



HNa 2 ^> 

ks 



*x(k) (4 



Bks C Aks - C Aks CB\cs 



(4) 



In the absence of interlayer Coulomb interaction and any 
phonon contribution, the total (volume averaged) c-axis 
current density would be equal to cr*(a;)£ ,ext (w) with 
a*(u) = (dbi/d)abi(w) and E cxt the external field, d is 
the c-axis lattice parameter. 

(ii) Response of a c-axis polarized infrared active 
phonon mode. In the absence of any electronic 

contribution, the current density can be expressed as 
u*{uj)E c ^{lo) with 



a*(u) = <7 p h(u;) = -iueox 



tdia-ph + dAjj (uj) 
i(u + iS) 



(5) 



where x is the phonon polarizability, Kdia-ph the term 
describing the diamagnetic contribution, and Ajj (w) the 
phonon current-current correlator. The term Kdia-ph is 
equal to — Eouip, where up is the phonon plasma fre- 
quency, up = l/(a 2 cfe ) (EpM^ 
e ^j ( e )p,z and m M are the effective charge, the c-axis com- 
ponent of the polarization vector and the mass of the /U-th 
ion. The correlator is given by 



Here 



Ajj(u) 



iNa 2 



dt< 



'<[&(*)'#h(°)]W*)- (6) 



wo the phonon frequency, a^ and a are the phonon cre- 
ation and annihilation operators, respectively. The cor- 
relator can be further expressed in terms of the retarded 
phonon propagator D(ui): 
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(8) 



In the noninteracting case, 

D{u) = D (u) = 2u a /(u 2 ~u 2 + iuS). 

(hi) Case of the electronic contribution coexisting with 
the phonon one. The total current density is given by 
a*{u)E eyA with 



a*(u) = 



Mia-bl 



+ Kdia-ph +dEjj((j) 



i(uj + iS) 



(9) 



where K dia _ b i = (du/d) Kdia-bh 



-jjM 



iNa 2 



die 1 ' 



x <[J b P ,W 
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^i(O) 



-ip p h (o)])^W (io) 



and J?, — {dbi/d)j?y In the absence of interlayer 
Coulomb interaction and electron-phonon interaction, 
the mixed correlators (involving both J^ and j? h ) 
would vanish and we would obtain simply o~*(u>) — 
(du/d)a-u(v) + cr ph (o;) with cr b i from Eq. (fT|) and er ph 
from Eq. §5§. The interactions will be shown to lead 
both to a renormalization of Cbi and cr p h and to nonzero 
values of the mixed correlators. 

(iv) Formulas used to describe the final state interac- 
tions. Let us begin our considerations with the simple 
scheme of Fig. 1 (b) of CBM. The difference between 
the local conductivities Obi and <7i n t leads to a differ- 
ence between the current densities jbi and jmt result- 
ing in a charge redistribution between the CuC>2 planes, 
that can be characterized by the excess surface charge 
density p of the bottom plane (the density of the upper 
plane is —p). The electrostatic energy due to p is given 
by Na 2 dbip 2 /(2eo) and the corresponding contribution 
to the Hamiltonian by Eq. (25) of CBM, i.e., 



Na 2 d b i „ 2 



(11) 



where p is the density operator, that can be expressed in 
terms of the quasiparticle operators as 



P 



2Na- 



Yl ( C Aks C Bks + 4ks C Ak S ) • (12) 



Here j^ h is the phonon paramagnetic current-density op- 
erator, 



7 P 

Jph 



i 1 Uujo ^— v 



.(e) 



ft z 



(a* — a) 



(7) 



Next we derive the formula for the coupling term iJ c _ p 
describing the interaction between the surface charge 
densities of the planes p and —p and a phonon. We begin 
with the textbook formula for the energy U of a polarized 
medium in a homogeneous electric field: U — —V ■ P ■ E, 



where V is the volume of the medium, P the polarization 
and E the field. In the present case, V is the volume of 
N unit cells, V — > Na 2 d, P is the polarization due to the 
phonon, P — > P, 



P 



1 
'aid' 



V 2w ^^ x/rru. 



(13) 



and E is the electric field due to p and — p, E —> p/eq. 
In addition, the product has to be multiplied by a fac- 
tor denoted by £, representing the degree of overlap of 
E and the microscopic dipoles {E is nonzero only in the 
bilayer unit); £ = 1(0) for phonons associated with vi- 
brations inside the bilayer unit (outside the bilayer unit) 
and (« 1/2 for phonons involving mainly the oxygens of 
the Cu02 planes. By combining the expressions for V, 
P, and E we obtain 



H„ 



Na 2 d - 
-e Pp- 



(14) 



(v) Modification of the response by Hq and H c - p . The 



function Ejj defined by Eq. (flOl) can be written as the 
sum of four contributions: xa [the term involving J£(t) 
and j* (0)], xb [&(*) and Jph(0)], Xc [J&(t) and &(«)], 

and xd b'phW ano - ^bi(O)]- These correlators can be 
expressed, using perturbation theory, in terms of those 
referring to the case of He = and H e - p = 0: II, j, 
IIjp , IIpj, Upp, Aj,-, Ajp, Apj and App. Two of the lat- 
ter, Ujj and Ajj, have been defined above, see Eqs. §S§ 
and ([5]), the others are defined similarly, for example 
A jP M = (tJVoV^n, d«e-*(^ h (f),P(0)])e(t). Fig- 
ure 1 shows the Feynman diagrams representing the cor- 
relators [in a) and b)] and the interaction vertices [in c)]. 
Also shown are the factors to be assigned to the ver- 
tices when calculating the contributions of the diagrams. 
Figure 2 shows the series of the diagrams corresponding 
to xa , Xc , Xc and Xd ■ The series can be summed up by 
standard techniques^ and the results are 
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i-ip npJ ' 
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(15) 



where II' , is the expression corresponding to the re- 
peating segment of the diagrams shown in Fig. 1 (d) , 



n' 



pp/pj 



^0 \«bl£o 



pp/pj ■ 



(16) 
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FIG. 1: Feynman diagrams representing the correlators de- 
fined in the text [a) and b)], the interaction vertices corre- 
sponding to He of Eq. (Til) and H c - P of Eq. fTJ} [c)], and 
two segments that appear in the series shown in Fig. [2] [d)]. 



(vi) Simplifications resulting from the gauge invari- 
ance. The sum in the numerator on the right hand side 
of Eq. © containing ftdia-blj ^dia-ph and the four corre- 
lators of Eq. (|T5|) can be considerably simplified by using 
the relations: 



IF, 



A 



iluso 



"Xbl 



Pi 



<bi 



iloeq 



IIjp , Upp = — Xbl 



-x 



Aj-p , A 



'bl 



so 



pp 



7 X ' 



(17) 



that follow from the requirement of the charge conser- 
vation formulated in different gauges of the electromag- 
netic potentials or equivalently from the requirement of 
the gauge invarianc o 26 i 35 . We obtain 

ws ... (dbi/d)xu + X + (1 - 2QxbiX MR x 
cr (w) = -ia;e - . (18) 

i - Xbi[(o/dbi)rx - 1] 

(vii) Tota/ electric field, from cr*(w) £o c(w). The to- 
tal electric field P consists of the external field E cxt and 
the induced field E lnd due to the charge fluctuations de- 
scribed by p and the phonon polarization described by 
P, E ind = (d hl /d)(p/e ) - (P/e ) = j/{iuJSo). We have 
used the continuity equation, icup = j^i and the rela- 
tion between P and the total phonon current density j p , 
—iwP = jp, j is the total current density. The c-axis 
conductivity a(cj) defined by a(uj) = j(uj)/E(uj) is equal 
toa*(uj)[E CVLt {uj)/E(uj)},i.e., 



a(u) 



er*(u/) 



1 - io-*/(ujs ) 



-IL0E{)- 



(du/d)xbi + X + (1 - 2£)xwx 



[1 - (du/d)]xu - X ~ [1 " 2£ + (d u /d)e]XuX 

(19) 




XD = 



FIG. 2: Series of the diagrams corresponding to the quantities xa, Xb, Xc, and xd defined in the text. The sum of the 
contributions of the diagrams is to be multiplied by d^/d 2 in case of xa and by dbi/d in case of xc and xd- 



(viii) Considerations of screening by high frequency 
processes lead to a slightly modified formula: 



o-(w) =r.h.s.of(19) 
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Xbl 


x^ — 
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(20) 

where £oo is the interband dielectric function. The total 
dielectric function is given by e(u>) — e^ + ia(co)/(ojso). 
It is an easy exercise to show that the final formula 
[ (fT9|) or (|20"|) ] is precisely the same as the one result- 
ing from the single phonon version of the phenomeno- 
logical multilayer model. 36 . The important point is that 
is has been derived here on microscopic grounds, using 
perturbation theory, with the electron-phonon coupling 
constant determined by simple electrostatics. Within the 
phenomenological model, the response of a phonon to the 
applied field is governed by the eigenvector-dependent lo- 
cal electric field. Here, the eigenvector enters the equa- 
tions through the coupling constant. Our result adds 



credibility to the central claim of Ref. |lfj, that at T° ns 
the spectral weight of the low-energy component of the 
(local) intrabilayer conductivity <7bi begins to increase, 
possibly a manifestation of precursor superconductivity. 
Future studies involving the present formalism in con- 
junction with the local conductivities resulting from mi- 
croscopic calculations such as those reported in CBM 
(rooted in the Fermi liquid theory) or those of the kind of 
Ref. [33 based on a phenomenological ansatz^ related to 
the resonating valence bond picture are needed to achieve 
a further progress. 
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